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Abstract 

A misfit parameter is used to characterize the degree of frustration of or- 
dered and disordered systems. It measures the increase of the ground-state 
energy due to frustration in comparison with that of a relevant reference state. 
The misfit parameter is calculated for various spin-glass models. It allows one 
to compare these models with each other. The extension of this concept to 
other combinatorial optimization problems with frustration, e.g. p-state Potts 
glasses, graph-partitioning problems and coloring problems is given. 

PACS numbers: 05.50.+q, 75.10.Nr, 75.50.Kj, 75.50.Lk 
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It is well established that two ingredients are necessary to characterize a spin glass: frus- 
tration (competition among the different interactions acting on a certain magnetic moment) 
and disorder, see for reviews. But up to now a quantitative description of frustration 
seems to be incomplete. Toulouse |2| has introduced the function $ = ]l(c) hj which mea- 
sures the frustration effect in a local region of a lattice, where c indicates a closed contour 
along the Ijj = +/ or —I bonds. However, this function cannot be simply generalized to 
other spin-glass models, especially, it is not suited to models without underlying lattices. 
Frustration has an effect on ground-state energy and entropy. This can be easily seen start- 
ing from a (unfrustrated) ferromagnetic system by replacing +1 bonds by —I bonds with 
increasing concentration p ||. The ground-state energy increases up to a critical concentra- 
tion p c . Near p c the ground-state entropy starts to increase. This is reflecting the fact that 
the problem to find the ground state becomes a problem of combinatorical optimization with 
a large number of optimal and nearly optimal solutions. The aim of this letter is to use just 
the energy increase due to frustration as its global measure. 

Firstly this concept was used to characterize the frustration effect in an amorphous Ising 
model with antiferromagnetic short-range interactions [|J. A misfit parameter 

\E id \ - \En\ 



m 



(1) 



\E U \ 

was introduced, where E is the ground-state energy of the frustrated system and E ld is the 
ground-state energy of a relevant unfrustrated reference system. The latter can be obtained 
by replacing all negative bonds by positive ones. For the ±7 spin glass the relation to 
Toulouse's frustration function can be seen by the expression given by Barahona || 

Eo = — \hj\ + 2 \hj\ , 

< ij > unsatisfied (2) 

edges 

where the first term represents E td and the second one the numerator of (1) having in mind 
that the restricted sum over unsatisfied edges is correlated to Toulouse's function by the 
total string length at minimal matching of elementary plaquettes with $ = — 1. 



The misfit parameter m of Eq. (1) is used to characterize the frustration in Ising zig-zag 
chains in dependence on the chain length || and in a neural network model [/J. It has been 
generalized to quantum systems || |J and to define the local misfit and the misfit of a 
cluster of spins ||. 

Now we introduce a modified misfit parameter. For a given state i of a system it is defined 
by 

fii = n(Ei) = J _ ™ , (3) 

max min 

where Ei is the energy of the state i. E^ in and E^ ax describe the minimal and maximal 
ideally possible energy values, respectively, where 'ideal' refers to the assumption that all 
local energies yield a minimal (maximal) contribution to the total energy. For spin glasses 
these energies have to be calculated assuming that all bonds are satisfied (nonsatisfied). 
Although, in general, E^ in and E l ^ ax do not represent necessarily energies of a real system, 
they often can be identified with energies of a special reference system. In any case, they 
represent lower and upper bounds for the possible energy range of the considered frustrated 
system. Therefore \Xi is restricted to the interval between and 1. We define the misfit 
parameter for a system by the misfit of its ground state /io = I^(Eq). 

To clarify the term 'ideal' energy, we will discuss the misfit parameter of a spin glass. As 
mentioned above, the minimal and maximal ideal energies correspond to a Active state, 
where all interactions are satisfied and nonsatisfied, respectively: 

E id . =-E id =-V (A) 

^mva ^max / , ryl > V*/ 

<ij> 

where the sum goes over all interactions. Obviously, E^ in is the ground state energy of an 
unfrustrated reference system, which can be obtained by replacing all Uj by their absolute 
values. The misfit parameter can be calculated from Eqs. (3) and (4) as 

2 ( J2<ij> \ ) 

Compared to Eq. (1) we get m = 2/xo- For the ±J spin glass /io is the fraction of nonsatisfied 
bonds in the ground state [[U| . 



For the well-known Sherrington-Kirkpatrick (SK) model |TTJ the minimal ideal energy be- 



longs to the ground-state energy of a reference system, in which the probability distribution 



of the interactions is a Gaussian one folded about zero 112 1 In this case Eq. (4 ) gives 



EZJb = -K a Jb = -V2/7T « -0.798, (6) 

where b denotes the total number of bonds in the system. Eq. (6) leads to /j,q = for the 
mean-field solution in JO]]. In other words, due to the mean- field approximation the frustra- 
tion in the system vanishes and the resulting system is a Mattis-like spin glass. The misfit 



value fj,Q for Parisi's improved replica solution JT3[ is given in Table I together with a collec 



tion of data for various spin-glass models and related combinatorial optimization problems. 



Derrida |14| has considered the random-energy model as an approximation to spin-glass 
models and has calculated lower bounds for the ground-state energies in any dimension. For 
the ±7 spin glass on a d- dimensional hypercubic lattice this approximation yields reasonable 



misfit values (/io > 0.11 and 0.17 for d = 2 and 3, respectively, and /xq > 0.5 — 1/y 2d/ln2 in 
the high-dimension limit d — » oo, which is lower than the lower limit of the fully frustrated 
±1 system (see Table I)). Otherwise, for the symmetric Gaussian model on a square lattice 
the ground-state energies of the random-energy model are lower than E^ in . 
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TABLES 

TABLE I. Misfit parameter fj,Q for various spin-glass models and other combinatorical opti- 
mization problems in the ground state. 
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d = 2,3,4 








0.3125 


6d hypercubic 


exact 


m 




0.5 - l/(2Vd) 


d > 8 hypercubic 
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0.5 - 1/V2d7r 


d > 8 hypercubic 
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0.417 


6d fee 


exact 






0.5 - l/(2d) 


d-dim. fee 


exact 


m 


p-partitioning (±7) 




Bethe lattice 


analytical approach 


m 


p = 2 b 


0.074 


z = 3 








0.230 


z = 8 








0.257 


z= 10 






p = 3 


0.105 


2 = 3 








0.317 


z = 8 








0.354 


z = 10 






p-coloring (±7) 




Bethe lattice 


analytical approach 


m 


p = 3 





z = 3 








0.032 


z = 8 








0.058 


z = 10 






p-state ±1 Potts glass 




Bethe lattice 






p = 2 


0.080 


z = 3 


MC and annealing 






0.134 


z = 4 








0.1975 


z = 6 








0.236 


z = 8 








0.265 


z = 10 






p = 3 


0.0068 


z = 3 


analytical approach 






0.172 


z = 8 








0.204 


z = 10 







a extrapolated from 1/N scaling 

b The misfit for the corresponding (p = 2)-coloring problem results in the same /xq values. 
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From Table I the various effects of dimension, coordination number, distribution and 
range of interactions and number of states per spin variable on the frustration can be seen 
in a quantitative manner. 

At least for small coordination numbers and dimensions their increases result in additional 
constraints and therefore in increasing Lio values. The effect of different coordination num- 
bers can be seen both for the two-dimensional ±7 Ising spin glasses with different lattice 
structures and for the p-state ±7 Potts glass, the p-partitioning and the p-coloring problem 
with different z values. The influence of different dimensions can be studied by comparing the 
results for two and three dimensions. The comparison between two- and three-dimensional 
±7 spin glasses and spin glasses on a Bethe lattice with the same number of nearest neigh- 
bors z shows that liq is stronger influenced by the coordination number than by the spatial 
structure and dimension. However, it can be seen also that llq is lower in the Bethe lattice 
than in higher netted lattices. 



Analytical expressions for Lio(d) are given for hypercubic fully frustrated systems 0] and d- 
dimensional antiferromagnets with triangular plaquettes at least in the high-dimension 
limit. For finite d the results for fully frustrated systems are proved as upper bounds for 
systems with equal probability of +7 and —7 bonds. 

In systems with a Gaussian distribution of interactions the energy can be decreased by 
choosing and frustrating that bond with the lowest strength in a plaquette. Therefore the 
misfit for such systems is smaller than for comparable systems with a ±7 distribution. 
Other relations between the parameters of a model and the resulting frustration can be inves- 
tigated by using the p-state Potts glass, the p-partitioning and the p-coloring problem. With 
an increasing number of colors p for the nearest neighbours of a site in a p-coloring problem 
the chance increases to give all neighbours another colors and therefore the frustration in the 
system decreases. On the other hand, if the number p of subsets in the p-partitioning prob- 
lem increases, the problem becomes more complicated and restricted, leading to an increase 
in frustration. In Potts glasses the number of possibilities to avoid frustration increases with 



rising number p of states per spin and consequently liq decreases. As outlined in |29| the 



p-state Potts glass can be understood approximately as an intermediate system between the 
p-partitioning and the p-coloring problem, which correspond to a ferromagnetic and an an- 
tiferromagnetic Potts glass with special magnetization constraints. Therefore, for the same 
p the fj,Q values in the Potts glass are smaller than those of the partitioning problems but 
larger than those of the coloring problems. Small deviations for the case p = 2 are due to 
different methods. 

Summarizing the results, we have introduced a global misfit parameter generalizing the frac- 
tion of unsatisfied bonds in the ground state of the short-range ±7 spin glass to other spin 
glass models and to other systems with frustration. It measures the influence of frustration 
on the ground-state energy and, for the first time, it allows to compare various spin glass 
models quantitatively. It can be applied to systems without an underlying lattice structure. 
An advantage over the former parameter (1) is its invariance against any linear scaling of 
energies. It means, e.g., that additional self energy terms leave the misfit unchanged. 
The presented concept can be applied to other systems with frustration at least for such 
problems, for which the cost functions can be transformed linearly without additional con- 
straints. This can be understood in terms of the existing transformations between various 
problems using the NP-completeness. But often, e.g. for the traveling salesman problem, 
such a transformation is accompanied with new global constraints. We will focus on this 
topic in a forthcoming paper. 

The presented parameter only refers to the energetic aspect of frustration. It is an open ques- 
tion whether a similar parameter can be found for the entropic characterization of frustrated 



systems. A preliminary answer for the ±J models is given by Vogel et al. ||20|| . These authors 
have calculated the fraction of bonds, which are satisfied in all ground states. The difference 
between unity and this fraction can be used as a global entropic measure for frustration. 
Generalizations are under consideration. 

We benefit from discussions with A.R. Ferchmin, A. Hartwig, K.-H. Hoffmann, A. Mobius, 
P. Polaszek, H. Rieger, E. E. Vogel and J. WeiBbarth. This work is supported by the DFG 
(project no. Ko 1416). 
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